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Superconductivity in the Hubbard model on a square lattice near half filling is studied using 
an optimization (or correlated) variational Monte Carlo method. Second-order processes of the 
strong-coupling expansion are considered in the wave functions beyond the Gutzwiller factor. Su- 
perconductivity of d^2_j,2-wave is widely stable, and exhibits a crossover around U = Uco ~ 12t 
from a BCS type to a new type. For U ^ Uco iU ^ U^o), the energy gain in the superconducting 
state is derived from the kinetic (potential) energy. Condensation energy is large and cx exp(— t/J) 
[tiny] on the strong [weak] coupling side of Uco- Cuprates belong to the strong-coupling regime. 
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PACS numbers: 74.20.Mn, 71.10.Fd, 71.30.-|-h 

In a superconducting (SC) transition, conventional 
BCS superconductors follow a low-frequency sum rule of 
optical conductivity cri(a;) 0. However, recent experi- 
ments have shown that cuprate superconductors violate 
this sum rule . This violation implies a gain in kinetic 
energy {K) in the transition, because the sum of oiIlo) 
is proportional to —K Such kinetic-energy-driven 
superconductivity sharply contrasts with that of the 
conventional BCS superconductors, where the transition 
is induced by the lowering of potential energy jsl . 

The Hubbard model on a square lattice, H = Ht + 
'^u = E(.,j> (cLcj> + c]^Ci^a'j+UJ2j nj]'nji, is of- 
ten used as a simple model which probably seizes the 
essence of cuprates 6] . In spite of its importance, reli- 
able knowledge is limited particularly in the intermedi- 
ate and strong coupling regimes. It is still controversial 
whether superconductivity is realized in this model. In 
the strong coupling region, the Hubbard model is mapped 
to t-J-type models, where the d2.2_j,2-wave superconduc- 
tivity is concluded by exact diagonalization |7| and vari- 
ation methods . On the other hand, in the weak- 
couplin g re gion, many quantum Monte Carlo (QMC) 
studies [lOj came to negative conclusions for U/t — 2- 
4. Unfortunately, QMC is ineffective in larger-?/ regimes 
due to the negative sign problem. In contrast, RPA cal- 
culations fluctuation exchange approximations Il2ll . 
renormalization-group [isl and perturbative |l5L llGj] 
studies concluded d2.2_j^2-wave superconductivity. Be- 
sides, variational Monte Carlo (VMC) studies using the 
Gutzwiller projection argued that an antiferromagnetic 
(AF) order prevails widely near half filling, and narrows 
the SC region 

The purpose of this letter is to resolve the above dis- 
crepancy, and to explore the possibility of the kinetic- 
energy-driven superconductivity in the two-dimensional 
Hubbard model. By carefully studying the wave func- 
tions with vital improvement on the Gutzwiller projec- 



tion, it is found that the d-wave SC state is stabilized 
even in the weak coupling region, but its energy gain is 
too small to be observed in QMC. Furthermore, we find 
a crossover at Uco ~ 12t, over which the SC transition is 
induced by the lowering of kinetic energy. This indicates 
that the high-Tc superconductivity should be understood 
in the context of str ong correlation. 

A VMC method hj is useful to our purpose for its 
applicability to any U/t. Generally, a Jastrow-type func- 
tion, 5* — V^, is used as a trial state, where <I> signi- 
fies a one-body (Hartree-Fock) state, and V a correla- 
tion factor. For V, the Gutzwiller (onsite) factor [20j . 
Vg — rijll ^ i^^ 9)^j1^jl]^ been often chosen for its 
siniplicitv. Although Vg is successful in t-J-type models 
1^ |2j] , it brings about unfavorable results for the Hub- 
bard model Il9[. which were expected to be remedied by 
improving Vn [2 21 123] . Since the relationship between the 
Hubbard and t-J-type models is, Ht-j 
[2^. which yields 
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one can make improvements by applying the strong cou- 
pling expansion, e"''^, to "^/g {— "Pq^)- 
Along this line, we employ "S^q — 



Vq^ with < /i < 1 in this 
is an asymmetric projection 
less-than-half filling (n < 1, n: 



)rk [23. 



Q3 
for 



work 123. Here, 

electron density), where 

'^j — = (1 ~ '^iT)(l ~ ^^'^ ''' runs over all 

the nearest-neighbor sites. This factor takes account of 
virtual states in the second order of the strong-coupling 
expansion, as explained in Fig.^ As a one-body function 
<&, we use the following: Fermi sea $f for a normal state, 
a mean- field solution <I>af(A) |2^ for an AF state, and a 
fixed-n BCS state with d3:2_y2-wave symmetry $sc(^i C) 
for a SC state Here, A and C are variational param- 
eters, although they denote gaps and chemical potential. 



respectively, in the mean-field approximations. 

Since our trial functions have up to four parameters 
to be op timized, the use of optimization VMC proce- 
dure [28| is practically indispensable. Special treatment 
for the present calculations is that after the convergence 
of optimization, we continue iteration for a while, and 
average the data obtained during this excess process. 
Thereby, the accuracy in energy is markedly increased, 
sometimes to the order of 10^^ t. We used lattices of 
L X L sites (L — 6-16) with closed shells and periodic- 
antiperiodic boundary conditions, and collected typically 
2 X 10^-10^ samples, keeping the acceptance ratio 0.5. 
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FIG. 1: Weight assignment of Vq. A closed circle indicates 
a doubly-occupied (d-) site, (a) A configuration appearing 
for U/t = 00. (b) A d-site arises by a single hopping process 
from (a) (a virtual state in the second order), (c) A d-site 
exists, and two hopping steps are needed to reach from (a) (a 
virtual state of the fourth order). 

First, we compare the total energy E/t between ^I/q 
and *G [FigEl^a) for n = 0.88]. In each state, E/t is 
noticeably improved by the correlation factor Pq, par- 
ticularly in large-?7 regimes. In Fig[21[b), the magnitude 
of the improvement is depicted for some values of n. The 
SC and normal states have large values which increase 
as n approaches 1, whereas the value of the AF state is 
small and decreases as n — > 1. The improvements in the 
former states are caused by the lowering of kinetic en- 
ergy Et = {Tit) (not shown), because Vq promotes the 
hopping between neighboring sites by making the con- 
figurations like Fig^b) advantageous. In contrast, in 
the AF state, the binding effect between a d-site and an 
empty (e-) site is already included in 'I'af; the number of 
broken antiparallel-spin bonds increases in a Neel back- 
ground when an e-site goes away from its partner d-site. 
Thus, the effect of Vq on <i>AF is small and reduced as 
n ^ 1. This result shows that the correlation factor Vq 
is essential for the normal and SC states. Henceforth, we 
focus on vpg. 

Next, we consider the energy difference between the 
normal and the ordered states (condensation energy), 
AE = - £;°'-dcr^ .^j^i^,;^ plotted in FigEI At 

half filling, the AF state is dominant for any value of 
U as anticipated, and has a peak at J7 = ^ 8t, 

which originates in a metal-insulator (Mott) transition 
of the normal state (2^. Phase-transition behavior was 
found in ^'g^, although does not give the lowest 
variational energy for n = 1. As shown in Fig|3| AE for 
the SC state abruptly increases at C/ = Usi ~ 6-5t, and 
the increase there becomes sharper as L increases (not 
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FIG. 2: (a) Comparison of total energy among various func- 
tions. For instance, 'Q/normal' indicates 'I' — Vq^f- (b) 
Energy improvement by on ^I/g in each state for several 
values of n. Solid (open, half-solid) symbols indicate the d- 
wave (AF, normal) state. For n ^ 0.85, the AF phase is not 
stabilized for any U/t. 
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FIG. 3: Energy difference (condensation energy) between the 
d-wave (AF) and the normal states, estimated by ^'q. Solid 
(open) symbols denote the SC (AF) state. Additional lines 
for n = 1 and n = 0.88 are oc t/U (dashed), and exp( — (7/t) 
(dash-dotted), respectively. Shown in inset is magnification of 
the smaIl-[/ region for n — 0.88. Additional lines are oc (U /t)"^ 
(dashed) and exp(— t/i7) (dash-dotted). An arrow for n = 1 
indicates Usi/t, and those for n < 1 U^o/t at which l\E has 
a maximum. 



shown). Note that this behavior is inherent in the SC 
state, and iJ"™™ is smooth near U^i. Around t/si, the 
quasiparticle renormalization factor Z for ^E'q^ vanishes 
(FiglBJ, and the small-g behavior of the charge structure 
factor N{q) changes from linear to g^-like (not shown). 
Therefore, we consider that a superconductor-insulator 
transition takes place, and ^'g-" becomes insulating for 
U > Usi- Actually, the parameter g exhibits a cusp, and 
A and fi drastically change at U = Usi- The insulating 
behavior with AE oc ^U is consistent with the aspect 
of t- J-type models H ll IH^ . 

When n deviates from half filling, the behavior of AE 
is no longer a transition, but becomes a crossover. The 
inset of FigO shows a close-up of AE of the small-C/ re- 
gion for n = 0.88. It is found that AE is very small 
for U ^ 6M (- IQ-H at U = 4t), and behaves mildly 
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like a power-law function [31| . It is probable that due to 
this imperceptibly weak superconductivity, many QMC 
studies for U/t = 2-4 [l^ overlooked signs of SC order. 
At [/ 6.5t, an exponential-like rapid increase occurs 
which again originates in , and then AE reaches a 
maximum a,i U — Uco as indicated by arrows in Fig|21 
[3^ . As we will explain shortly, the SC properties for 
U > Uco qualitatively differ from those for U < Uco- For 
U > Uco, AS is nicely fitted by the form exp{~U/t) [or 
exp(— J)], which means that the effective attractive in- 
teraction is J = At^ /U; a viewpoint from i- J-type models 
is justified. In this region, the SC state is fairly stable 
to considerably large U/t. As n decreases, the difference 
between the two regimes tends to be vague. 

Before the SC properties, let us discuss the phase dia- 
gram. As shown in the inset of Fig|3l due to the sudden 
increase of AE^^ for n < 1, the stable phase switches 
from SC to AF, but the SC phase is retrieved at a larger 
value of U (FiglSJ. The obtained phase diagram in the 
U/t-n plane is shown in Fig^fa). The AF state is stabi- 
lized near half filling where the nesting condition is satis- 
fied, and occupies the maximum range of n at C/ ~ 6t. If 
we use instead of the AF region becomes larger. 
This implies that the improvement by e^*'^ is more effec- 
tive on the SC state than on the AF state. We expect 
that the higher-order improvements enlarge the SC re- 
gion. 
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FIG. 4: (a) Phase diagram within 'i'Q near half filling. The 
boundary for U/t ^ 6 is well fitted by 1 — n oc (U/tf, which 
looks similar to the result of a renormalization-group theory 
[1^ . For systems with n < 0.84, the d wave is always stable 
for U/t > 0, as illustrated by gray arrows, (b) Momentum 
distribution n(fc) of the d-wave state. M point is indicated by 
an arrow. 

We turn to the SC properties in Figure [S] shows 

the difference of kinetic energy AEt {— and 
potential energy AEij. Apparently, the mechanism of 
energy gain changes at ?7 ~ Uco- For U ^ Uco, the energy 
gain is due to the lowering of Eu {AEu > 0), as in the 
BCS superconductors. Inversely, for U ^ Uco, kinetic- 
energy gain occurs {Et is lower in the SC state) , agreeing 
with cuprates. This qualitative change is a crossover and 



not a phase transition. These aspects do not depend on 
n qualitatively, although the magnitude of both AEt and 
AEjj decreases, as n decreases. The lowering of Et has 
been found also by a dynamical cluster approximation 
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FIG. 5: Difference of kinetic and potential energies between 
the normal and SC states for several values of n. The absolute 
values tend to increase, as L increases (not shown). 
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FIG. 6: Quasiparticle renormalization factor in the nodal 
direction for the normal (open) and the d-wave (solid) states. 
The inset shows the difference of Z between the normal and 
the d-wave states. 



The momentum distribution function n{k) gives an- 
other evidence for the crossover. Since there are nodes in 
the d3.2_j^2-wave state, n{k) has a discontinuity at A; = fcp 
in the y = ±x (F-M) direction [Fig^b)]. The magni- 
tude of this jump represents the renormalization factor 
Z, which is related to the inverse of the electron effective 
mass. The estimated Z for 4-"°™ and ^-^^ are plotted in 
Fig|Sl For U ^ Uco, the normal state has a larger value 
than the SC state, whereas this relation is reversed for 
U ^ Uco (see also inset). For U > Uco, electrons are hard 
to move (or heavy) in the normal state, hampered by the 
strong correlation. On the other hand, in the SC state. 
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the coherence somewhat reheves the suppressed raobihty, 
and thus the energy gain occurs in Et- 

This kind of crossover from weak- to strong-coupHng 
regimes is not restricted to the present case. The at- 
tractive Hubbard model exhibits similar behavior in the 
s-wave SC transition for any electron density Hence, 
the kinetic-energy-driven superconductivity arises, not 
because the model has repulsive interaction, but the in- 
teraction is strong enough. The VMC calculations for the 
t-J model showed that the exchange (hopping) energy is 
lowered (raised) in the d2,2_j,2-wave SC state (See Table 
I in ref.Q, and also [l^l)- Since the kinetic part of the 
Hubbard model induces the exchange term (J = At^ /U), 
the results of i-J-type models are compatible with the 
present results in the Hubbard model for U ^ Uco- 

The magnitude of condensation energy in our calcula- 
tion is larger than that obtained experimentally [s^ • We 
anticipate that further refinement on the wave functions, 
especially on vl/norm^ leads to a better outcome. 

In summary, we have studied the two-dimensional 
Hubbard model near half filling, based on evolved VMC 
calculations. It is shown that the da;2_j,2-wave SC state is 
stabilized in the Hubbard model if the variational states 
are improved. Furthermore, we find a crossover at around 
Uco/t = 11-13. For U ^ t/co, the SC state is similar to 
the conventional BCS state, but for U ^ Uco, the SC 
transition is caused by lowering the kinetic energy. Re- 
cent experiments of ui (w) for cuprates 2] accord with the 
characteristics of this strong-coupling region. Extended 
description will be given elsewhere |l8l |. 
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